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Abstract — Multicast transmission has several distinctive traits 
as opposed to more commonly studied unicast networks. Spe- 
cially, these include (i) identical packets must be delivered suc- 
cessfully to several nodes, (ii) outage could simultaneously happen 
at different receivers, and (iii) the multicast rate is dominated by 
the receiver with the weakest link in order to minimize outage 
and retransmission. To capture these key traits, we utilize a 
Poisson cluster process consisting of a distinct Poisson point 
process (PPP) for the transmitters and receivers, and then define 
the multicast transmission capacity (MTC) as the maximum 
achievable multicast rate times the number of multicast clusters 
per unit volume, accounting for outages and retransmissions. Our 
main result shows that if r transmission attempts are allowed in 
a multicast cluster, the MTC is 6 [pk^ log(fc)) where p and x are 
functions of r depending on the network size and density, and k 
is the average number of the intended receivers in a cluster. We 
also show that an appropriate number of retransmissions can 
significantly enhance the MTC. 

I. Introduction 

Multicast refers to the scenario whereby a transmitter needs 
to send a packet to multiple receivers. In a wireless network, 
this creates a two-edged sword. On one hand, the broadcast 
nature of wireless transmission assists multicast; but roughly 
uncorrelated outage probabilities at each receiver (due to spa- 
tially distinct fading and interference) require retransmissions 
that cause interference and waste. Multicast is an important 
aspect of sensor and tactical networks, and increasingly in 
commercial networks where streaming is supported. However, 
the literature on multicast is minuscule compared to unicast. In 
this work we attempt to investigate the fundamental throughput 
limits of multicast transmission and we develop a metric 
based on spatial outage capacity which we term multicast 
transmission capacity (MTC). 

In order to characterize the MTC in a wireless network we 
propose a multicast network model in which each transmitter 
has an intended multicast region (called a cluster) where all the 
intended receivers are uniformly and independently scattered, 
and hence a Poisson cluster process can be reasonably used 
to model the transmit-receiver location statistics. The active 
transmitters are modeled as a stationary Poisson point process 
(PPP) and their associated receiver nodes in the cluster are 
also a stationary PPP, as shown in Fig. [T] In other words, 
each cluster is randomly located in the network and comprises 
a multicast session. This paper will develop interference and 
outage expressions for this model, and analyze some important 
cases of the network model and design space. 
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Fig. 1. Multicast transmission model in a 2-dimensional ad hoc network: 
Each transmitter (triangle) and its intended receivers (small circles) are 
indicated by the same color in a cluster of radius s. 



Some previous works, such as Pl-fS), have made significant 
progress in studying the multicast or broadcast capacity. For 
example, in [1] the protocol model is used where source nodes 
and their multicast destinations are randomly chosen. The 
multicast capacity is defined as the sum rate of all multicast 
flows and it is obtained as a function of the number of 
multicast sources. In IS), the multicast capacity under the 
protocol model is defined as the transmission rate summed 
over all of the multicast traffic flows in the network. Its scaling 
characterization is obtained by the number of receivers in each 
multicast session. In |7|, the physical model and a stationary 
PPP of the nodes in the network are considered. It showed that 
the broadcast capacity is a constant factor of the computed 
upper bound when the number of nodes goes to infinity under 
a constant node intensity. 

The multicast and broadcast capacities in the prior works are 
not investigated from the multi-receiver outage point of view 
and thus their scaling results cannot provide us retransmission 
guidelines for outage reduction and capacity enhancement. The 
multicast capacity problem here is studied from a multicast 
outage perspective. The main issues we would like to clarify 
are that how many multicast sessions which interfere each 
other can coexist if all receivers in a multicast session need 
to receive the packet from their transmitter, and when retrans- 
mission is beneficial or detrimental to the multicast capacity. 

In this paper, we introduce the MTC in a d-dimensional 
network, which is defined as the maximum achievable multi- 
cast rate per transmission attempt times the maximum number 



of the coexisting clusters in the network per unit volume 
subject to decoding delay and multicast outage constraints. 
The idea of defining MTC is from the transmission capacity 
framework originated in lH. The decoding delay constraint 
here means a transmitter can multicast a packet to all of its 
intended receivers up to t & N+ transmission attempts and 
multicast outage happens when any of the intended receivers 
in a cluster does not receive the information multicasted by 
their transmitter during t attempts. 

Main Contributions. In this work our first contribution 
is introducing a first way to define the MTC with multicast 
outage, and propose some new cluster-based definitions for the 
largeness and denseness of a network in order to characterize 
the scaling behaviors of the MTC under different conditions. 
We show that the MTC scaling can be expressed in a general 
form of 8(pfc^ log(/c)) where p and x are given in TableUfor 
different network conditions, k is the average number of the 
intended receivers in a cluster. From the scaling characteriza- 
tions, we know retransmissions have a significant effect on the 
MTC and certain number of retransmissions could enhance it. 
In addition, we also show that the MTC scaling is not affected 
by the fading model of communication channels. 

TABLE I 

Main Results on Multicast Transmission Capacity 
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II. Network Model and Preliminaries 

The multicast network model should be with tractability 
of finding the necessary information for MTC. The following 
Poisson cluster model is proposed because it not only captures 
the multicast behavior in the network but also makes all 
derivations easier. 

A. Multicast Transmission Model 

In the network, each transmitter has a multicast cluster 
of equal volume and its receive nodes in the cluster suffer 
aggregate interference from a Poisson field of transmitters. 
Specifically, we assume that the network is operating a slotted 
ALOHA protocol and the distribution of the transmitting nodes 
in the network is a stationary Poisson point process (PPP) $t 
of intensity At. As shown in Fig. [T] any transmitter Xi e $t 
has its own intended multicast cluster TZi where all of its 
intended receivers are uniformly and independently distributed 
and they also form a stationary PPP of intensity Ar. 
Note that each cluster could contain other transmitters and 
unintended receivers except its own transmitter and intended 
receivers. 

Accordingly, the multicast transmission sessions in the net- 
work are essentially the Poisson cluster process Zi = UXi, 
i.e., each transmitter is a parent node associated with a 
cluster of receive daughter nodes. The cluster processes {Zi] 



corresponding to different transmitters {Xi] are assumed to be 
independent so that the superposition of all clusters yields the 
resulting cluster process $ — [Jx £<i>t intensity X — kXt 
where k — is the average number of the intended 

receivers in each cluster assuming all {Tli,yi G N} have the 
same volume The distribution of the intended receiver 
nodes in each cluster is assumed as a marked PPP denoted by 
= {(Ytj,Htj) : Y^j £ TZi,j G N}, where Hij is the fading 
channel gain between transmitter Xi and its intended receiver 
Yij. Similarly, the distribution of the transmitters in the net- 
work is also a marked PPP, i.e., $f = {(X,, {-ffy }), i,j G N} 
where Hij denotes the fading channel gain between transmitter 
Xi and the receiver Ybj located in cluster 7?.o. All the fading 
channel gains are i.i.d. with probability density function (PDF) 
fH{h). 

Without loss of generality, the MTC can be evaluated in 
the reference cluster 7?.o whose transmitter Xq is located at 
the origin. We condition on this typical transmitter Xq result- 
ing in what is known the Palm distribution for transmitting 
nodes in the d-dimensional Euclidean space f9l. It follows by 
Slivnyak's theorem [91 that this conditional distribution also 
corresponds to a homogenous PPP with the same intensity 
and an additional point at the origin. The signal propagation 
in space is assumed to undergo path loss and fading. The path 
loss model between two nodes X and Y used in this paper is 



\X-Y\ 



\X^Y\ 



if > 1 

otherwise, 



(1) 



where |X — y | denotes the Euclidean distance between nodes 
X and Y, and a > d is the path loss exponent. The reason 
of using the model in ([TJ is because the model | • does 
not behave well in the near field of each transmitter and it 
thus leads to an unbounded mean of the shot noise process. 
The Nakagami-m fading modeQ is adopted here because 
it covers several different fading models, such as Rayleigh 
fading, Rician fading and no fading, etc. By using this model, 
we can observe if different fading models affect MTC or not. 

The intended multicast region TZi of transmitter Xi is 
confined by B{Xi,s) which is a d-dimensional ball centered 
at transmitter Xi with radius s > 1. Each receiver in is 
able to successfully receive the message if its SIR is greater or 
equal to the target threshold /3. That is, receiver node Yj G $r 
is "connected" to the typical transmitter if 



>/3, 



(2) 



where all the transmitters are assumed to use the same transmit 
power, the network is interference-limited, and It is the 

' is the Lebesgue measure of TZi . For example, if TZi is a rf-dimensional 
ball of radius s then pLr = /^u s**, where /lu = \/7r^/r(l + d/2) is the 
Lebesgue measure of a rf-dimensional unit ball |'9 |. 

^In this paper, the Nakagami-m fading channels are always assumed to 
have unit mean, unit variance and m S N4.. 

^Here means ■ The subsubscript "0" of <^ro is dropped for notation 
simpHfication. Since all of the following analysis is based on the nodes in the 
reference cluster 7?-o, the subscript or subsubscript of some variables will 
not be explicitly indicated if there is no ambiguity. For instance, Yj and Hj 
in TZq actually stand for Y^j and Hoj, respectively. 



aggregate interference at receive node Yj and a sum over the 
marked point processes. Namely, 

It^ H,,\\X,~Y,\\-^, (3) 

Xie*t\{^o} 

which is a Poisson shot noise process, and Hij is the fading 
channel gain from transmitter Xi to receiver Yj in TZq. Since 
$t is stationary, according to Slivnyak's theorem the statistics 
of signal reception seen by receiver Yj is the same as that 
seen by any other receivers in the same cluster. Thus It can 
be evaluated at the origin, i.e., ^ can be rewritten as It = 

J2xie<S>t\{Xa} ^iW-^iW 

Suppose the decoding delay is up to the lapse of r transmis- 
sion attempts for a transmitter. The connected receiver process 
for the i-th transmission is denoted by 

1.,, = e : H,^ > mrit} , (4) 

where {Hj.} are i.i.d. for all i g [I,-- - ,t]. Also, let $c 
be the connected receiver process at the rth attempt, i.e., it 
is the set of all intended receivers in a cluster connected by 
their transmitter during the decoding delay, and thus it can 
be written as $c — Ur=i In other words, the connected 
receiver process can be described by a filtration procesfl 

B. Multicast Transmission Outage 

The transmission capacity of an ad hoc network introduced 
in |[8l is defined based on point-to-point transmission with an 
outage probability constraint e E (0, 1), and is given by 

Ce = 6At(l-e), (5) 

where b is the constant transmission rate a communication link 
can support (for example, about log2(l + /?)), and Af is the 
maximum contention intensity subject to an outage probability 
target e. However, ^ cannot be directly applied to multicast 
because the multicast rate will be affected by /i^ and Xt, and 
the outage of a multicast transmission is not point-to-point but 
point-to-multipoint. 

Since no desired receiver can be assumed to be dispensable, 
a reasonable way to define outage is when any of the intended 
receivers of a transmitter does not receive a multicasted packet 
during a period of time up to the decoding delay. Thus, 
a multicast outage event of each multicast cluster can be 
described as £ = {$c C The probability of 8 can be 
characterized by the intensity of the connected receivers during 
the lapse of r attempts as follows: 

p[5] ^ 1 - P[{$A *c} = 0] 

= l-exp|-^ (A, -Ac(r,T))M(dy)|, (6) 

where is a d-dimensional Lebesgue measure. Since all of the 
intended receivers are uniformly distributed in 7?.o, (01 can be 
rewritten as 

P[£] = 1 - exp {-^ir {Xr - Kr[Xc{R, t)])} < e, (7) 

*A filtration process means C ' ' ■ C , and for any set A. C 
TZo, 3>cj (^) — > "3?r(^) almost surely as j — oo where denotes the 

random number of point process <1> enclosed in set A. 



where R e [0, s] is a random variable whose PDF is 
//j(r) — d ^-p-- The outage probability in (|7]i cannot exceed 
its designated upper bound e which is usually a small value. 

Remark I: When a transmitter has a fixed number of the 
intended receivers in a cluster, the probability of multicast 
outage based on our definition is too complex to be calculated 
due to the spacial correlation of interference ifTOl . Since $r is 
a stationary PPP and $c is a nonhomogeneous PPP (this point 
will be proved in Section HIIl l. we can obtain the result in (|6]l 
so that the multicast outage probability can be easily carried 
out via finding E/j [Ac(i?, r)] in (|7]l. Since Ei^[Xc{R,t)] is a 
monotonically decreasing function of Af, At can be reached 
by the lower bound on Efj[Xc{R, r)] obtained from (|7]l- 

C. Definitions of MTC, Largeness and Denseness of Networks 

The MTC defined in the following is similar to the idea of 
the transmission capacity defined in Q, but there are some 
subtle differences. Since the multicast outage probability is 
upper bounded by a small e, a Taylor expansion gives Af (e) — 
Ae(e) + O(e^), where A^ is a reasonable approximation of Af 
for small e. For simplicity, we will focus the analysis on A^. 
In addition, the transmission rate 6 in (|5]l becomes a multicast 
rate, which is not necessary equal to a constant. 

Definition 1: The multicast transmission capacity with the 
multicast outage probability defined in (|6]l for small e is 
defined as 

^-6A,(l-e), (8) 

T 

where h is the maximum achievable multicast rate on average 
for every cluster. 

Cc essentially gives area spectral efficiency of cluster- 
based multicast transmission. The following definitions of 
largeness and denseness of a network will be needed to acquire 
the scaling characterizations of the MTCs in the subsequent 
analysis. They are defined based on the large average number 
of the intended receivers fc in a cluster 

Definition 2: (a) We say a network is "large" if the volume 
jir of a cluster in the network is sufficiently large such that for 
fixed Xr we have fc 1. (b) If the intended receiver intensity 
is sufficiently large such that for fixed volume we have 
fc ^ 1, then such a network is called "dense", (c) A "large 
dense" network, it means that clusters in a network have a 
sufficiently large size as well as receiver intensity; namely, 
Xr (X fj,r and thus fc ^ 1. 

III. The Receiver-Connected Point Process 

During the allowed r transmission attempts, the intended 
receivers in TZq connected by transmitter Xo form a receiver- 
connected point process whose intensity is the necessary infor- 
mation to estimate the multicast outage probability. Since $c 
is a filtration process and upper bounded by $r (as explained 
in Section Ill-Al l, the connected receiver intensity in TZq is an 
increasing function of r as shown in the following lemma. 

Lemma 1: Consider the stationary PPP in the reference 
cluster TZq. If a transmitter is allowed to transmit a packet 
up to T times and all channels are Nakagami-m fading, then 



<i>c is a nonhomogeneous thinning PPP and its intensity at 
r € [1, s] is lower bounded as follows. 



1 - 



1 - 



(/3r")^'(™-i)(/3r") | , (9) 

A. A.(,,oo)/^)^ (10) 



where Ai is defined in the following lemma for moment 
generating functional of stationary independent PPPs. 

Proof: See Part lAl in the Appendix. □ 
Lemma 2: Let = {{X^^ , H,^ ) ; e 6(0, r) nlR'', r > 
1 , j e N} be a stationary marked PPP of intensity Ai for all 
i e [1,2,- •• and {Hi } are i.i.d. Nakagami-m random 
variables with unit mean and variance. Suppose $i has a 
Poisson shot generating function 



is defined as /; 



E 



which 



where a > d. If 



are independent and ^ = ^, then the sum of the Poisson 
shot generating functions, i.e., / = Ei=i ^i' has the following 
moment generating functional for A — Ei=i ^ 1^^++ 

and (/)2 € (0, m r"): 



>C/(0i) 

X/(02) 

where 

Ai((Ai,r) : 

A2(02,r) : 



E [e-*i^] = exp Ai (01 , r) A) , (1 1) 
E e^^^l =exp(M„A2((/.2,r)A), (12) 



m/01 (1+t)"- 
m/<f>2 



dt, 



di. 



(1-t)' 

Proof: Omitted due to space. See lITTI if interested. □ 

IV. Multicast Transmission Capacity 

In this section, we study the MTC when transmitters are 
multicasting to all of their intended receivers in a single- 
hop fashion. Using the receiver-connected intensity found in 
Section ini] and the multicast outage probability defined in (|7J, 
the maximum contention intensity is found as follows. 

Theorem 1: Suppose dTji is upper bounded by e and the 
maximum decoding delay is t transmission attempts. If fc > 
-^rrc, then the maximum contention intensity is 



A. 



e 



(13) 



where p = ^ e{T + 1). 

Proof: See Part iBl in the Appendix. □ 
Remark 2: The scaling function &{■) of A^ in (fT3] l only 
contains the "controllable" network parameters such as iJr, 
\r, e and /?, which means their values are adjustable if 
needed. The parameter m of Nakagami fading, is a channel 
characteristic and usually uncontrollable and thus not involved 
in the scaling result in Theorem [T] Hence, the scaling of the 
maximum contention intensity is unaltered no matter whether 
the channels in the network are fading or not. 



If a unicast planar network without retransmission is con- 
sidered (i.e., d = 2 and k = t — 1), in ( fTsT l will reduce 
to the previous results discovered, i.e., Ae = 8 p2/a ^ ■ 
In 1 8 1, for example, the maximum contention intensities of 



FH-CDMA and DS-CDMA are 9 



and e 



respectively, where M is the channel number of FH-CDMA 
and the spreading factor of DS-CDMA. It is easy to check 
that these two results coincide with ours here by considering 
for FH-CDMA and 4 for DS-CDMA. In addition, the 



M 



longest transmission distance in a cluster is s and we know 
Ae = Q{s~^) and so is the network capacity, which also 
coincides with the results in JSl IIT2I . 

Considering a capacity-approaching code is used, the max- 
imum achievable multicast rate b that is acceptable for all 
intended receivers is the following Ergodic channel capacity 
evaluated at the boundary of a cluster: 



6 = E 



log 1 



It 



(14) 



where H„ 



The bounds on b are given 
in the following lemma. 

Lemma 3: There exists a 5 G (0, 1) such that the bounds 
on the multicast rate b in ( fT4l ) can be given by 

(51og(l + -V ) <^<logl 1 + ) +0(1). (15) 



\ Xt J \ ^J^rXt , 

Proof: Omitted due to space. See ifTTI if interested. □ 
Remark 3: The bounds on b in ( fTSl l are not affected by 
channel fading since the fading effect has been averaged out, 
and they indicate b is significantly dominated by the aggregate 
interference from the transmitters in the cluster 

Scaling Law of MTC. According to Theorem[T]and Lemma 
[3] we found that the multicast rate b is 0(log(fc)/T) for any 
network conditions if A^ is achieved and k is sufficiently large. 
The MTCs in a network without receiver cooperation can be 
concluded as follows, (i) For a dense network, Ae = 8 
since fir is fixed and /c ^ 1. By the MTC definition, we 



know Cf 



(ii) If the network is large, then 



A, = 8 (fc^S^) and thus C, = 
large dense network, A^ — 8(/ij.^ 



^ \ fc(l+l/^) 



(iii) For a 

and A; 1. So A^ is 
e (^pr^ k^i^)] and thus = (^pk'i^) log{k)j . In 
summary, the MTC here can be expressed in a general form 
as follows: 

= 8(pfc-log(fc)), (16) 
where x has been given in Table |T] 
Appendix 

A. Proof of Lemma Q] 

The Laplace functional of the stationary PPP $ for a 
nonnegative function 5 : M'' — > IR+ is defined and shown 
as follows ifTsI : 



E 



exp 



- 4<i giX) -S-idX) 



3-»W) A.Md^)] 



Since the Laplace functional completely characterizes the 
distribution of the point process, we can find the inten- 
sity of by looking for £$^(5). Recall that $c = 
{Yj e <^r ■■ max,g[i^....^] Hj. > /3||Y,||"/f} and {iJ^J are 
i.i.d. \/i E [1,2,- ■■ ,r]. Let 1^(2^) be an indicator function 
which is equal to 1 if a; € A and 0, otherwise. The Laplace 
functional of $c for g{Y) = g{Y)l^^{Y) is given by ( fTTI i. 



00 , j 

i=0 



n 



l-E[U^iY,)] 



fi{dY,)---pi{dY,), 



■ exp 



(17) 



Also, for all 7' e [1, s] we know 



max iJj, > /3||F|r/t 



(*) 



where (★) follows from the fact that the temporal correlation 
of interference can be neglected for small A* fTQ\. So we have 



^<E>c (5) = exp - 



where Ac(r,T) = A,(l-(E/[^^H(/3r" (/^ + iVo)|/t)])"). From 
the above result we know that $c C $r is a nonhomogeneous 
PPP because its intensity Xc{r,T) is the intensity A^ of $r 
scaled by (fTST l. Since the channel gain is a Nakagami-m 
random variable, we know that 



E, [F^(/3r"/t)|/t)] 



> 1 - 



r(m, f3r°'uj) 



/7,(w)dw, 



/o r(m) 

where r(m,a;) = J^f^-^e-*dt and r(m) = r(m,0) = 
(m — 1)!. Also, if fw{w) is a probability density function of 
random variable W then we can have the following result. 

r{m,aw)fw{w)dw^{^ar^-^ ■ (18) 



According to Lemma|2] we can have Ci^{4)) with Ai(0, 00). 
Thus, 



E/ [FH{lir''It\It))] > 1 



r(m) d^"" 



-1 /£7j0)\ 



1 - (/3r")*("-i)(^r°). (19) 



Substituting ( fT9] l into Ac(r, r), then (|9]l can be arrived. 

B. Proof of Theorem Q] 

According to the outage probability (|7]i upper bounded by 



e, we know 



Efl[Ae(i?,T)] > A, 



log(l - e) 

Mr 



A.(i4)+e(.^), 



for sufficiently small e. Calculating the lower bound of 
Efl [Ac(i?,T)] from (|9]l and then using the Holder inequahty, 
we can show that 



Efl[(/3i?")*('"-i)(/3i?")] 
Note that 



> 1 



(20) 



G (0, 1) almost surely 
and thus its average approaches to unity when e/k is 
sufficiently small such that \f^Jk < e. That means 
(/? i?")^^™^^'' (/3 i?") is nearly equal to one almost surely 
and thus At = 0(e)- If k is sufficiently large, then we have 
exp(— /i„ AtAi) = 1 — At Ai + 9(e^). Substituting this 
expression into (/3 (^i?") can be reduced 

to 1 - /i„ At Ai nr=~i (1 - e/j) + 0(e')- Define Ai{^) ^ 

[Ai((/),oo) - Ai(0,aB)]nr=i'(l-e/j). Choosing ag < 1, 
we have 



E, 



(1- (;3i?")^'("-i)(/3i?"))^ 



< 



Another upper bound for the above equation can be obtained 
by (|7|i and (|9]l, and it should coincide with the above upper 
bound when At is maximized as shown in the following. 



A, = 



(21) 



flrf3i V^Ai(/?)' 

So ( fTSl l is obtained and the proof is complete. 
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